Exact solution of the ion-laser interaction in all regimes 
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We show that in the trapped ion-laser interaction all the regimes may be considered analytically. 
We may solve not only for different laser intensities, but also away from resonance and from the 
Lamb-Dicke regime. It is found a dispersive Hamiltonian for the high intensity regime, that, being 
diagonal, its evolution operator may be easily calculated. 
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I. INTRODUCTION 

The ion-laser interaction may be easily solved in the 
low intensity regime (LIR) [l|-l3], but besides the condi- 
tion that the laser intensity is much lower than the vibra- 
tional frequency, we set the condition that the detuning 
between the laser and the atomic transition frequency is 
an integer multiple of the vibrational frequency. Then 
some questions arise: Is it possible not to consider inte- 
ger multiples of the vibrational frequency? Is it possible 
to solve for high and middle intensities? 

Indeed, Moya-Cessa et al. [a| have shown that it is 
possible to find solutions for any set of parameters, i.e. 
in all the regimes. However the solutions are not general 
because the set of eigenstates found can not expand all 
possible (general) states. 

It has been shown already that for low intensities it 
is possible also to consider the ion micromotion [3|, and 
by using Ermakov-Lewis invariant methods 10| it was 
possible to linearize the ion-laser Hamiltonian when the 
micromotion was included [11]. Here we would like to 
show how it is possible to solve the interaction in different 
regimes, including high intensity and medium intensity, 
the method allows also not to consider multiple integers 
of the vibrational frequency. 



II. ION-LASER INTERACTION 

The Hamiltonian for the ion-laser dipole interaction, 
with no approximations can be written as 

H = i'h + -;^az + i^{fJ+ + a-) 



where v is the harmonic trapping frequency, LUa is the 
atomic transition frequency, uj]^ is the field frequency, 51 
the (real) Rabi frequency of the ion-laser coupling and 
77 the Lamb-Dicke parameter. The operators a^ and a 
are the creation and annihilation operators for the vibra- 
tional motion of the ion, and the ct's are the Pauli spin 
operators. 

By doing the transformation Tlip) with T__ = 



exp(- 



:uia+S 



Gzi) and performing the optical RWA [1^ we 



arrive at the well-known Hamiltonian 



Hion ^vn+ -a^ -I- il [cF+b{iiq) + a-D^ir]) 



(2) 



where D{ia) = e'"(''+° ) is the displacement operator, 
and (5 = cjf, — ujl the laser-ion detuning. 



A. Low intensity regime 

The low intensity regime is the well-known regime, 
where several effects like multi-phonon transitions, 
Jaynes-Cummings (JC) and anti-JC interactions may be 
obtained. To solve this regime, we follow first the tra- 
ditional approach. We start by using Baker-Hausdorff 
formula [15.] to factor the displacement operators in equa- 
tion (2) into a product of exponentials and consider 5 — 
kv, i.e. an integer multiple of v with /c = 0, ±1, ±2, ..., 
we then obtain 



Hion = vn+—crz+ne 



-.V2 (^ 



cr+e 



^-irja' ^-i-qa 



(3) 

Now we expand the exponentials of the annihilation and 
creation operators in Taylor series and get rid off the 
free Hamiltonians via a transformation to the interaction 




P(t) 




20 30 



FIG. 1: Plot of Pe{t) as a function of t for it = 0, i^ = 1, ^ = 
0.2 and rj = 0.1. Solid line represents the numerical (exact) 
solution, dashed line the solution from Section II and the 
dot-dashed line the solution for the dispersive Hamiltonian of 
Section III. 



picture to obtain the Hamiltonian 



(4) 



^-^ n\m\ ] 
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We use the fact that are in the LIR, v ^ Q. and make 
the RWA, i.e. we only keep time independent terms in 
the above Hamiltonian to end up with 



i/j = fie-'''/2 atfc(_j^)fe^_ 



"■ if(r;^) + ff.c 



(5) 
with L^ ' [x) the associated Laguerre polynomials of or- 
der (operator) h = a^ a. The Hamiltonian above is now 
readily solvable, so that we may find easily the evolution 
operator associated to it. 



III. OTHER REGIMES 

Although the atom-field and ion-laser interactions ap- 
pear to be physically and mathematically quite distinct, 
they are in fact exactly equivalent. The easiest way to 
see this is by using the transformation 



R 



6 2 £4 



(<T+-o-_)g-i^(o-|-at)((j^) 



FIG. 2: Plot of Pe(t) as a function of i for it = 0, fi = 1, 
!/ = 0.2 and 77 = 0.1. Solid line represents the numerical 
(exact) solution, dashed line the solution for the dispersive 
Hamiltonian of Section III . 



such that Hion = RB^on 



n 



-VLc 



R) 



)(^ 



{^^ 



(6) 



(7) 
Therefore we have linearized the ion-laser interaction 

in an exact way. In the following we will neglect the term 

2 
■^ because it only represents a constant shift of all the 

eigenenergies. 



A. Medium intensity regime (MIR) 

We now consider the case where the vibrational fre- 
quency is of the order of (twice) the field intensity (Rabi 
frequency). We also consider the Lamb-Dicke regime, 
i.e. 77 <C 1. For simplicity we will set 5 = to show 
the different possibilities we have now. However it is not 
difficult to produce effective Hamiltonians also in the off- 
resonance case. In this case the Hamiltonian ([7]) may be 
casted into 

'HMiR. = vn + Q.<7z + — {&+a + +a)a-) (8) 

which is a Hamiltonian that has been extensively studied 

MM- 



B. Low and high intensity regimes (HIR) 

We have shown in Section II how to solve for the LIR 
case. Here we will show a different method that is also 






valid for the HIR. 

By transforming the Hamihonian ([7|) with the unitary 
operators 



C/i 



^1 (a ' (7+ —aa- ) 



C/2 



o52(aa-+-ata-_) 



(9) 



with ^1,^2 ^ 1; we can remain up to first order in the 
expansion e^^^Be-^^ = B + ^[A, B] + ^[A, [A, B]] + ... « 
B + £^[A, B], so we obtain the effective Hamihonian 



n 



eff 



va)a + Vl<Tz - XionC^zia^a + -) 



K 

2' 



ttcTz ia) 



(10) 



We have used 



rjv 



2{v + 2VL) 



6 = 



rjv 



2{2Vl-v)' 



(11) 



We can see that in fact ^i , ^2 ^ 1 either in the LIR 
(in this case we have also to consider 77 <IC 1) or in the 
HIR (no constrain for rf) , which justifies completely the 
approximation for the above Hamihonian. For the reso- 
nant case, A: = 0, it becomes diagonal and we can solve it 
in an easy way. In Fig. 1 we show a plot for the probabil- 
ity to find the ion in its excited state. The three curves 
in the figure correspond to the exact case (solid line) , the 
solution form Hamihonian ([S]) (dashed line) and the so- 
lution for the dispersive Hamihonian pUj) . We can see 
excellent agreement among the three plots for the LIR. 
Now, for the HIR we show a plot in Fig. 2 for the numer- 
ical solution (solid line) and our solution from this Sec- 
tion (dashed line). Again it may be noticed an excellent 
agreement between both curves. We should stress that 
there is no other analytical solution to compare with, as 
ours is the first analytical solution in this regime (also in 
the medium regime). 

The new interaction constants in the effective Hamil- 
tonian ([TU)) have the form 



4f]2 _ 1,2 



S-qv"^ 



4f]2 



(12) 



In the resonant case and high intensity regime, fi ^> i/, 
it is easy to show that 



Xion ^ X 



high 



2?72;/2 1 
4" 1-^ 



2n ' 



(13) 



while in the low intensity regime, 51 <C i', we will have 
the same Hamihonian but x will change to 



Xlo 



-2rfVt- 



1 



4r!2 



-2rfn 



(14) 



If in Eq. (TTOl) we take the detuning 5 different from 
zero, we could get the usual blue and red side-bands in- 
teractions (see for instance Ref. |7l,i|). This is done by 
choosing the value 5 — ±z/. The only case in which 
we can obtain such regimes is the low intensity case, 
where one can perform the RWA to the Hamihonian (|10l) , 
which agrees with the usual procedure for obtaining such 
blue and red side-band regimes. The high intensity case, 
Q, ':^ V does not allow such side-bands because in the 
Hamihonian ()10p the interaction constants multiplying 
the different terms may be of the same order. 



IV. CONCLUSIONS 

We have shown that it is possible to solve analytically 
the ion-laser Hamihonian in different intensity regimes, 
from low to high. For the MIR we have casted the ion 
laser Hamihonian into a JCM Hamihonian (for the on- 
resonant case) that allows easy solution. For the HIR 
we have found a dispersive Hamihonian, which, being 
diagonal, it is direct to solve. We have found excellent 
agreement between the exact (numerical) solutions and 
our proposed solutions. 

We would like to thank CONACYT for support. 



[1] D.J. Wineland, J.J. Bollinger, W.M. Itano, F.L. Moore, 

and D.J. Heinzen, Phys. Rev. A 46 R6797 (1992). [8] 

[2] D. Leibfried, R. Blatt, C. Monroe y D. Wineland, Rev. 

of Mod. Phys. 75 281 (2003). [9] 

[3] J. Wineland et al., J. Res. Natl. Inst. Stand. Technol. 

103 259 (1998). [10] 

[4] D. Leibfried, D.M. Meekhof, B.E. King, C. Monroe, 

W.M. Itano, and D.J. Wineland, Phys. Rev. Lett. 77, 

4281 (1996). [11] 

[5] S. Wallentowitz and W. Vogel, Phys. Rev. A 59, 531 

(1999). [12] 

[6] M. Sasura and V. Buzek, J. 0/ Mod. Optics 49, 1593 

(2002). 
[7] H. Moya-Cessa and P. Tombesi, Phys. Rev. A61, 025401 [13] 



(2000). 

H. Moya-Cessa, D. Jonathan and P.L. Knight, J. of Mod. 

Optics 50, 265 (2003). 

W.P. Schleich, "Quantum Optics in Phase Space", 

(Wiley- VCH, 2001). 

M. Fernandez-Guasti, and H. Moya-Cessa, J. of Phys. 

A36, 2069 (2003); H. Moya-Cessa and M.F. Guasti, 

Phys. Lett. A311 1 (2003). 

J.M. Vargas-Martinez and H. Moya-Cessa, J. of Opt. B 

6, S618 (2004). 

We can perform the optical rotating wave approximation 

as the laser and atomic transition frequencies are much, 

much larger than the field intensity in experiments [j]. 

E.T. Jaynes and F.W. Cummings, Proc. IEEE 51, 89 



(1963). [15] W.H. Louisell, Quantum Statistical Properties of Radia- 

[14] B. W. Shore, and P. L. Knight, J. Mod. Opt. 40, 1195 tion, (Wiley, NewYork, 1973). 

(1993). 



